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Abstract 

We propose a model of a spatially modulated collective charge state (CCS) 
of superconducting cupratescJ. The regions of higher carrier density (stripes) 
are described in terms of one-dimensional (ID) interacting fermions and the 
regions of lower density as a two-dimensional (2D) interacting bosonic gas 
of dx2_y2 hole pairs. The interactions among the elementary excitations are 
repulsive and the transition to the superconducting state is driven by decay 
processes. Vibrations of the CCS and the lattice, although not participating 
directly in the binding mechanism, are fundamental for superconductivity. 
The superfluid density and the lattice have a strong tendency to modulation 
with wave- vectors (vr/a, 0) and (0, vr/a) implying a still unobserved dimerized 
stripe phase in cuprates. The phase diagram of the model has a crossover 
from ID to 2D behavior and a pseudogap region where the amplitude of the 
order parameters are finite but phase coherence is not established. We discuss 
the nature of the spin fluctuations and the unusual isotope effect within the 
model. 

PACS numbers:74.72-h, 74.80-g, 74.20.-z 
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I. INTRODUCTION 



The experimental evidence for charge and spin inhomogeneities in superconducting 
cuprates has been accumulating in the last few years!'!. It is believed that these inhomo- 
geneities are probably the result of strong competing electron-electron and electron-lattice 
interactions. The strong interactions can be hinted from the fact that cuprates are related 
to Mott insulators, not to metals. From the theoretical point of view the major problem has 
been the description of the doped Mott insulator. It is this state that finally evolves into a 
superconducting state. The origin of the inhomogeneities can be diverse. Disorder and the 
natural charge localization effects that occur close to an insulating Mott state are natural 
possibilities. Furthermore, there is an emerging point of view that inhomogeneities may be 
in the heart of the superconducting phase. The difficulty in describing the superconducting 
state may be related with the fact that an important part of the physics occurs in real space. 
This should be contrasted with a Fermi liquid description that is dominated by the Fermi 
surface in momentum space. If the inhomogeneities observed experimentally are essential 
for the description of the superconducting state then a dual description in terms of real and 
momentum space is certainly required. 

The search for a general principle or a special (maybe hidden) symmetry that facilitates 
the understanding of these materials is always desirable. It turns out that the intrinsic 
complexities of these systems are enormous. In fact, the experimental evidence seems to 
be that all the degrees of freedom, charge, spin, lattice, participate in a fundamental way 
in their physics. Unfortunately, we may have to sacrifice the simplicity of description that 
helped us to understand the basic physics of many systems in the past. Besides the strong 
local interactions that induce large charge fiuctuations, long-range interactions associated 
with the insulating state can suppress the same degrees of freedom. We believe that it is the 
interplay between these two forces that is fundamental for the understanding of the physics 
in these systems. In this paper we assume that the physics in cuprates can be divided at 
least into two length scales. On the one hand there is the short length (high energy) scale 
that is associated with the behavior of the electrons at atomic distances. This piece of 
physics is dominated by the strong local interactions that lead to the Mott insulator and 
it is well described by models of strongly correlated electrons such as the Hubbard or the 
t-J model. This is a "highly" quantum regime where quantum fiuctuations at atomic scales 
determine the physical response. Since the short length scale is usually smaller than the size 
of systems studied numerically, these studies can provide a good physical insight. Various 
analytical approaches to these strongly coupled models exist in 2DqD and although they 
can also provide insight into the main physical aspects of the problem they are unreliable 
in dealing with quantum fiuctuations. In order to acquire intuition one is forced to look at 
purely academic models that can be solved beyond mean-field. This is the case of the various 
studies of the Ising t — J model where quantum spin fiuctuations are suppressed by a large 
anisotropy in spin space!, the spin-fermion models where the spins are treated classically! or 
models with large number of components!. On the other hand, the long wavelength physics 
of the problem is dominated by long-range interactions such as the Coulomb, CasimirEl or 
entropic interactions that are essentially classical in originElS. Analytical treatments of 
these interactions have always been difficult even in statistical mechanical problems in ID. 
Therefore, from the theoretical point of view there are technical complications at all length 
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scales. Any theoretical treatment that intends to start from the microscopic picture from 
the beginning will encounter tremendous difficulties. 

In this paper we follow a semi-phenomenological approach. At the short length scales 
we have known for a long time that doped Mott insulators have a big "aversion" to charge 
homogeneitylii. This happens because there is a large loss of magnetic energy when charge 
delocalizes. Even before superconductivity was discovered in cuprate oxides the concept of 
strings and quantum confinement of holes was well understood: when a hole moves in an 
antiferromagnet it produces a string of overturned spins with an energy growing with its 
length0. When the hopping energy, t, is smaller than the antiferromagnetic exchange, J, 
strings are the dominant processes in antiferromaenets (although higher energy delocalizing 
processes called Trugman loops are also possibldl^) . When two holes are injected into an 
antiferromagnet there is a gain in magnetic and kinetic energy if they move together instead 
of separately. Strings are suppressed because the second hole heals the string generated 
by the first one. This picture of the diluted Mott insulator is generic. It can be obtained 
analytically in the Ising t — J model and is also observed in numerical simulations of the 
SU(2) models as welfi Because of the tendency to segregate chargej_Mott insulators are 
at the edge for phase separation into hole-rich and hole-poor regional^. Still to this date 
there is a heated debate about the border to phase separation in the phase diagram of 
the t — J mo More recently the concept of stripes in the t — J model has emerged 

due to density matrix renormalization group (DMRG) calculations by White and Scalapino 
(WS)El. These simulations show clear signs of charge order. In the WS picture when holes 
are doped into the antiferromagnet they first form pairs that condense into lines of charge 
with magnetic anti-phase domain walls (ADW). This condensation process, as shown in 
analytical approaches to the Ising t — J model£l and in other numerical studiesil, has to do 
with the gain in kinetic energy of a single hole due to the formation of an ADW (moreover, 
Trugman loops are supressed in this cas That is, for a finite linear density of holes 
(vanishing 2D density) the gain in kinetic energy by making a domain wall is enough to 
compensate for the loss of magnetic energy due to the formation of a magnetic defects. We 
should stress that we are not concerned with the problem of the phase diagram of the t — J 
model because this phase diagram most certainly cannot be the phase diagram of cuprate 
oxides. There are many important interactions in cuprates that are not included in the t — J 
model. We will argue below that coupling to the lattice, for instance, is relevant for the 
experimental phase diagram. We assume that the tendency of the diluted Mott insulator to 
form bound state pairs of holes is universal. 

Long-range interactions and their eventual screening are probably the key to understand 
the phase diagram of cuprate^ll. At low doping, because of the existence of a large charge 
gap in the Mott insulator, charge dynamics is suppressed and as a consequence dynamical 
screening as well. All the physics rests in the spin degrees of freedom. Thus, when a small 
concentration of hole is introduced into the Cu02 layers long-range interactions should play 
a major role. These long-range interactions will be effective in suppressing the tendency to 
phase separation as stressed by Emery and Kivelson (EK) and coUaboratorsil^ and will 
lead to the generation of a finite length scale, Ld-, associated with domain sizeii. In the 
absence of a lattice the long-range forces produce blobs or lakes of charge with characteristic 
size Lr,. It turns out, however, that in transition metal oxides the coupling to the lattice 
is strong and therefore the symmetry of rotations in real space is broken. Charge modu- 
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lation is therefore the final result of the tendency to phase separation and strong lattice 
coupling. The formation of this charge modulated phase is not the result of a Fermi surface 
singularity as in the case of ID charged density waves (CDW). In fact, the CCS state de- 
scribed here has strong similarities with the metallic 2D charge density wave states observed 
in dichalcogenides^ that have very little resemblance with the ID CDW states that are 
usually insulatingil. This state is probably better described as coming from the quantum 
melting of an anisotropic Wigner crystal (or Bragg glass)0 or an exotic CDW stateil. In 
classical statistical mechanics there are many analogues called liquid crystal phases: ne- 
matic, smectic, hexatic, etcil. These are phases with long-range orientational order in the 
absence of translational order. Another more mundane phase where translational symmetry 
is broken in 2D is an anisotropic membrane phase where translational order is broken along 
the principal vectors of the latticed. 

Let us consider a static charge modulated phase induced along the y direction in the 2D 
system with periodicity Na where a is the lattice spacing. The charge density, p(r), can be 
written as (we use units such that h = ks = I) 

p{r) = po + mr)) exp + c.c. (1.1) 

where po is the background density, and $ is the complex order parameter that can be 
rewritten in terms of an amplitude |$| and phase u: 

(<l>(r)) = |$|exp|-^|^n(r)| . (1.2) 

The lines of constant phase at wave number 271 /{N a) are describe by 

(l) = ^(y-u) = 2nn (1.3) 
Na 

where n = 0, ±1, ±2, .... In a inert background (such as a Fermi liquid) the local modulations 
of the density can lead to gaps at the Fermi surface but the system as a whole remains 
homogeneous. In a Mott insulator this is not necessarily so. In the regions of low charge 
density the Mott insulator is essentially untouched. In these regions it is energetically 
advantageous for the system to form bound states of holes with d^_j^2 symmetry (exactly 
like in a finite cluster or ladder) with a gain magnetic energ In the regions where 

charge density is large it is energetically more favorable for the system to have a gain 
of kinetic energyS that liberates the holes to move as single particles and create ADW. 
Thus, the formation of a charge ordered state has to be accompanied by a change in the 
spin structure of the system with the creation of incommensurate spin fiuctuations. The 
mechanism of gain of energy (kinetic or magnetic) depends strongly on the amount of charge 
density. Moreover, deformations of amplitude of the order parameter of the CCS are always 
energetically costly. Therefore, phase fiuctuations are the low energy excitations in such 
systems. For a classical smectic phase, for instance, the free energy is given by 

Fs = \jdh [B{d\\usf + K{dlusf 

= ^E^kKk|' (1.4) 
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where B is transverse stiffness of the smectic and K is the splay elastic constant. Here || 
and _L indicate the gradient of the function in the directions parallel and perpendicular to 
the ordering direction, respectively. Moreover, cj^k = -B/cj^ + Kk'^ is the classical smectic 
dispersion relation. A membrane phase, on the other hand, would be described by ujm,^ = 



-B|^A;j^ + Bj_k'j_ where B\\ {B±) are the longitudinal (transverse) compressibilities. In fact, 
because the charge modulated state is composed of electrons we expect this modulated phase 
to be quantum in nature, that is, a quantum liquid crystal. 

Notice that such a smooth charge distribution as proposed in (|1 . 1| ) is different from 
the stripe phases discussed in the context of insulating cuprates and nickelates. In in- 
sulating nickelates static charge order is observed in neutron scattering and many other 
experiment^. The key point is that higher harmonics of the fundamental Bragg peaks 
associated with static order are observed in these systems. It implies sharp, well defined 
and isolated domain wall structures in real space. In superconducting systems where in- 
elastic neutron scattering peaks are observed at equivalent positions there are no harmonics 
observed, even when the fundamental peaks are rather sharp. In Lai,48Ndo.4Sro.i2Cu04 
(LNSCO) where Tranquada and collaborators observed quasi-static peaks in neutron scat- 
tering there are no signs of higher harmonicJU. This experimental fact signals to a smooth 
variation of the charge-spin densities in the system. Thus, the idea of well-defined, non- 
interacting stripes is misleading in the context of superconductivity. Instead one should 
think of the "stripes" as a complex collective state that is driven by competing local inter- 
actions and cooperative long-range forces. 

The importance of the lattice degrees of freedom has been experimentally verified in es- 
sentially all cuprate superconductors and their insulating relativesll'll. At this point in time 
most of the theoretical approaches either focus on electron-electron interactions and overlook 
the importance of lattice degrees of freedom or mainly electron-phonon interactionJ^ dis- 
regarding the importance of the strong (short- and long-range) interactions in the problem. 
Some mean-field approaches, however, have stressed the importance of electron-lattice inter- 
action in the context of nickelatesiS. We believe that electron-electron and electron-lattice 
interactions are equally important because of charge neutrality. Charge neutrality implies 



that a charge modulated state such as the one defined in (|T]T|) has to be strongly coupled 
to the lattice. One would expect under general considerations that fluctuations of the CCS 
to appear in the phonon spectrum. Indeed, consider the classical elastic lattice free energy: 

^ k 

where wp^k is the lattice deviation from equilibrium position and cup^k the phonon dispersion 
relation. In the leading order the coupling between the smectic and lattice is quadratic in 
the displacements and can be written as 

Fc = J2 CkMkMp-k + c.c. (1.6) 
k 

where Ck are the coupling constants. The problem of the CCS plus phonons can be solved 
exactly by a simple linear combination of the displacement fields. There are two branches 
of excitations with frequency given by 
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leading to a splitting of the vibrational modes. Therefore, the phonon spectrum should be 
directly affected by the presence of a CCS. Thus, in dealing with the fluctuations of the 
CCS we have to consider the renormalizations of such fluctuations by the lattice. The full 
quantum mechanical problem can be quantized exactly like phonons in ordinary crystalline 
solids and the vibrations of the CCS can be described by 

= E^"AA>k. (1.8) 
k,a 

where ba^p {bl_^) is the annihilation (creation) operator for quantum vibrations of the CCS 
with momentum k in the branch a and energy fi^^k- 

Although the description of the collective state is rather simple because of its gaussian 
nature, the description of its internal degrees of freedom (associated with the short length 
scales) is more complex. In what follows we will make an ad hoc assumption that the internal 
charge excitations of this collective state can be divided into two main groups. Firstly, in 
the regions where the charge density is large (that is, given by the lines of constant phase 
in ( p.. 3D ) ADW exist due to the local gain in kinetic energy. These regions we call stripes. 
These high density regions are characterized by single particle excitations (not bound statesT) 
and they are essentially confined to ID lines because of the potential induced by stringsEj. 
This highly anisotropic electronic fluid should be interacting because of the phase space 
constraints and in the absence of tunneling between stripes it is described by a Luttinger 
liquid^. In the regions of low density (in the middle of the antiferromagnetic ladders) 
single particle excitations are suppressed and bound states of holes are energetically more 
favorable. The simplest of them is a bound state of a pair of holes with (ix2_y2 symmetryS. 
Larger bound state structures like quartets are unlikely to contribute because their quantum 
dynamics is exponentially suppressed (tunneling matrix elements decay exponentially with 
the number of particles). Moreover, we assume that this gas of bosons is essentially isotropic 
and only weakly coupled to the stripe fermions. The main reason of the weak coupling is 
related with the "string healing" process that generates the pairs in first place, that is, the 
bosons are essentially insensitive to the magnetic structure including the ADW. We argue, 
however, that it is the weak coupling between bosons and stripes that ultimately produces 
superconductivity in the system. 

From the electronic point of view the situation is illustrated in Fig.|T[ The Mott insulating 
state can be described by a filled lower Hubbard band (L.H.B.) and an empty upper Hubbard 
band (U.H.B.) that are separated by a large energy scale. In between these two bands there 
is a single electron band associated with the stripe fermions that we call the stripe band. The 
bound states of holes exist due to the transfer of spectral weight from the L.H.B. to a level 
above it with energy i?k/2 (where k is the momentum of the boson as we discuss below). The 
binding energy of the holes is the energy difference between the boson level and the top of 
the L.H.B.. The binding energy can be seen as mediated by the exchange of paramagnonJ^^I 
and therefore is of the order of the characteristic magnetic energy in the problem0'@. In the 
undoped system the characteristic energy is simply the exchange constant J but as doping 
increases the magnetic energy scales are reduced driving the system from the Mott insulator 
to a more ordinary Fermi liquid state. 
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The measure of the magnetic energy is the spin stiffness, ps, of the magnetic background. 
A finite spin stiffness produces the confining potential for the charge carriers of the form0: 

Vc{y) = ^\y\ (1.9) 

(modulus Na) as shown in Fig.0(a) (Vdy + Na) = Vc{y)). The total potential as seen by 
the holes is a superposition of the atomic potential of the lattice, Va(i/) iVA{y + a) = Va(|/)), 
and the magnetic confining potential of the strings. In Fig.^(b) we show in a simple Kronig- 
Penney picture the result of the superposition of these two potentials. When the temperature 
is larger than pg the holes are essentially deconfined. Since the hole-pairs only exist at 
temperatures below ps it is important to estimate this energy scale. In order to do so let us 
consider the situation in Fig.|^(a) where the exchange within the antiferromagnetic regions is 
J but across the stripe it is J' < J. At long wavelengths the problem maps into a spatially 
anisotropic Heisenberg model with exchanges Jx and Jy as shown in Fig;^(b)ii. While J'/ J 
is a short length scale property and depends on the microscopic detailsEJ the ratio 

« = ^ (1.10) 

Jy 

depends on the long wavelength properties and determines the region of stability for an- 
tiferromagnetic orderii. Moreover, the effective spin stiffness of the magnetic background 
depends directly on a. Assuming that Jy = J in Fig.^(b) (so that the exchange along the 
stripe is not modified by the presence of the stripes) it is easy to show thatil 

Ps{a) = JS^^/a. (1.11) 

This choice for the exchange constants was used in ref.@ and seems to explain well the data 
m LNSCO§. In the low hole doping regime (x < 0.02) another choice has to be made since 
the magnetism is isotropic^. It is clear from (|1 . 1 1|) that a decrease in a reduces the spin 
stiffness and can lead to a quantum phase transition to a spin liquid state@. The main 
problem is how to relate a of the effective model of Fig.|^(b) to the microscopic model of 
Fig.^(a). We use a simple scheme and compare the classical ground state energies of the 
two problems. In the case of Fig.^(a) the Hamiltonian is simply 

H = J2J^,A■^J (1.12) 

i,j 

where Jjj = J in the antiferromagnetic ladders and Jjj = J' across the ADW. Let us rewrite 
J' = •yJ where < 7 < 1 is a microscopic quantity that depends on details of the problem. 
With this parameterization it is clear that the classical ground state energy can be written 

as 

Eo = Eaf{J) + (1 - l)JS^N,trrpesNs (1.13) 

where 

Eaf{J) = -2JS^N,pins (1.14) 
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is the classical energy of an isotropic antiferromagnet with Ngpins spins and exchange J. Ns 
is the number of sites in each direction (iV^ is the total number sites), Ngtripes = Ns/N is 
the number of stripes in the system {N is the separation between stripes in lattice units), 
and Nspins = N"^ — x Ngtripes is the number of spins not residing in stripes [Nspins = 
Ar2(i _ 1/A^)). Thus, from ( [LT3D one gets. 



-2 + —-^. (1.15) 



A^2J52 N 

On the other hand, for the effective model shown on Fig.^(b) with the same number of sites 
we would have 

Eo = -{J, + Jy)S^Nl = - J(l + a)S^N^ (1.16) 
where we used (|1.1CI| ). Comparing (|1.15|) and ( p..l6|) we find 



a = l-^. (1.17) 



This result implies from (|1.11| ) that 



p.(iV) = JS^l-^. (1.18) 

Them the spin stiffness vanishes at a critical distance Nm between stripes given by 

Nm = 3-7. (1.19) 

Since < 7 < 1 we see that 2 < A^^ < 3 (quantum fluctuations increase the value of A^m^)- 
The binding energy of the hole-pairs is proportional to ps and no bosons can exist at zero 
temperature when < A^^. Moreover, we identify a temperature scale T*{N) ^ Ps{N) 
above which no bosons exist {T*{Nm) =0). So we relate the formation of the hole pairs 
as the "pseudogap" energy scale observed in many different experiment^. Our theory, 
however, breaks down when T > T*. 

In the non-interacting picture (with no hybridization between bosons and fermions) we 
can have two possibilities. In Fig.0 the boson state is empty of holes (filled with electrons) 
and the stripes are partially filled up to the chemical potential energy p. Because bosons 
and fermions are in thermodynamic equilibrium we work in the grand-canonical ensemble 
and keep the chemical potential fixed by letting the number of particles fluctuate. Consider 
the case when 



P > E^/2 . 

The boson state is unoccupied and there are no bosons in the system. As the chemical 
potential is reduced (more holes are introduced) the stripes empty and further reduction of 
the chemical potential pins the energy at the boson level since as more holes are added they 
can only produce bound states. Thus, there is a continuous transfer of spectral weight from 
the L.H.B. to the boson level. On the other hand, if the binding energy of the holes increases 
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nothing happens until fj, = Ek/2. At this point there would be bosons and stripes co-existing 
with each other. When the binding energy is increased further the chemical potential follows 
since bosons become converted into stripe fermions. Once again the chemical potential is 
pinned at the boson level. Finally when the boson level reaches the top of the stripe band 
the stripes are completely filled and the ADW disappears. In fact, because of the loss of 
kinetic energy the domain wall disappears even before this limit is reachedlli. Therefore, we 
assume that 

2/i - > . (1.20) 

Moreover, it is clear from this picture that the spin degrees of freedom that are responsible 
for the large magnetic response in these systems leave on the L.H.B. and are separated in 
energy from the charge degrees of freedom. Therefore, the magnetism can be effectively 
"traced out" of the problem since it will only lead to kinematic renormalizations of the 
various parameters (in other words, the spins follow the charge). 

The mechanism for superconductivity discussed in this work requires the coupling of 
stripe fermions via the exchange of bosons. One can think of this mechanism as the exchange 
of stripe "pieces". Coherence between the fermions requires the exchange of real bosons. 
The simplest mechanism for exchange is the decay of the bosons (since they are composite 
particles) into fermionic degrees of freedom at the stripes. This kind of mechanism can occur 
when two systems with very different ground states are separated by an interface. In fact, it 
was proposed long ago that a mechanism of this sort could generate superconductivity at a 
metal-semiconductor interface (also called exciton superconductivity)!^. Although this kind 
of proposal has generated controversy in the past^ there are good indications that they may 
be good candidates in the case of cuprate^^l. Moreover, the process described here is similar 
to the "proximity effect" mechanism proposed by EK if we neglect retardation effects due to 
the boson motionil^. One of our main results is: if the stripes static, that is, if we disregard 
the fluctuations of the CCS, this process is suppressed! The "proximity effect" mechanism 
cannot happen with static stripes and real bosons. Therefore for coherence to be attained 
we have to introduce fluctuations of the CCS. This makes our model radically different from 
the proximity effect proposed by EK since it involves distortions of the CCS and therefore 
of the lattice as well. Moreover, unlike BCS, the phonons are not part of the binding 
mechanism that is driven by the exchange of composite fermion pairs. Vibrations, however, 
are fundamental for the superconductivity. A natural consequence of the mechanism is that 
fluctuations of the order parameter are necessarily coupled to the fluctuations of the CCS 
and therefore vortices are coupled to dislocations of the stripe array. Because the superfluid 
density is low, the interactions between topological defects ultimately determines the phase 
diagram. Moreover, we claim that static stripes (as the ones observed obtained in Hartree- 
Fock solutions of the t — J and Hubbard model ME or DMRClll) should be insulating due 
to a ID CDW instability along the stripe direction. Diagonal stripes do not couple to the 
bosons because they are oriented along the nodes of the boson wave-functions. 

Among other things we explain why phonon anomalies that have been observed in neu- 
tron scattering! occur exactly at the same position in the Brillouin zone where angle re- 
solved photoemission (ARPES)i observes the opening of the pseudogap (that is, at (vr/a, 0) 
and (0,7r/a)). Moreover, we show that the lattice distortions associated with the phonon 
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anomalies are associated with the modulation of the superfluid density perpendicular to 
the stripes. The theory also predicts an isotope effect and the co-existence of commensu- 
rate and incommensurate spin fluctuations. The critical temperature for superconducting 
order, Tc, is determined by the interplay between topological defects associated with the 
superfluid bosons (vortices and anti- vortices) and distortions of the CCS (dislocation loops). 
The transition at finite temperatures is in the 3D XY universality class. We show that 
the amplitude of the order parameters is finite at temperatures above Tc up to T* in the 
pseudogap region. Moreover, at T = there is a quantum phase transition as a function 
of the separation between stripes, A^, that is directly related to the hole concentration x 
[N ^ l/(2x)). While the order parameters become finite for < Ngp long-range order 
is only attained at = A'^c < A^sp- This quantum phase transition is in the 2D XY in a 
magnetic field universality class. Thus, there is a crossover region at T = 0, A'c < A^ < Nsp, 
where topological defects prevent long-range order to develop. 

The paper is organized as follows: in the next section we introduce the bosonic bound 
state of holes; in section |T| we present the Luttinger liquid representation of the stripe 
fermions; in section |IV| we discuss the nature of the coupling between bosons and stripe 
fermions and argue that vibrational degrees of freedom should enter explicitly; in section 
|V] we solve the mean-field equations for the problem and determine the mean-field phase 
diagram; in section |V| we discuss the phase fluctuations of the superfluid-superconducting 
state and how they relate to the phase fluctuations of the CCS; in section [Vll| the nature of 
the high temperature transition is discussed; section [VIII] contains the results for the zero 
temperature transition; section 11^ contains our conclusions. 



II. ANTIFERROMAGNETIC BOSONS 

Consider the problem of two holes in an antiferromagnet. The energy is minimized by 
the formation a bound state that moves freely through the systemislil. One can define an 
operator that creates such a state at the top of the antiferromagnetic vacuum that is given 
by 

= 4^ II ^k(q)Ct,k/2-qCi,k/2+q (2.1) 
V J q 

where Co-,k destroys an electron (creates a hole) with momentum k and spin projection a (f 
or j) in the antiferromagnet. S = (A^^a)^ is the area of the system. The wavefunction of 
the pair, \E'k(q), is normalized {J2q |\E'k(q)P/5' = 1) and in principle depends on the relative 
and total momentum of the pair. This wave-function can be obtained variationallyii or by 
solving the Bethe-Salpeter equation for the binding of two holesil. The dependence of the 
wave-function on the total momentum is due to the fact that on the lattice the symmetry 
of the bound state varies with its center of mass momentum and depends strongly on the 
microscopic details. For the SU(2) t — J model it is known that the da;2^2-state has the 
lowest energy!^ while in the Ising t — J model the p-wave is the ground stateilii. The matter 
of the fact is that the s-wave state is always the state with highest energy and the reason 
is fairly simple to understand: the strong local repulsion requires the wave-function of the 
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pair to be centered at different sites. In other words, s-wave bound states are suppressed by 
the magnetism. Here we assume the pairs to have d^^-yi symmetry. 

It is easy to see that operator defined in is not completely bosonic because of its 
composite nature!^. Since the electrons obey anti-commutation relations {{ca,^, c^, ,^,} = 
6q^q'6cr,cr') it is easy to show that commutation relations of the bosons is given by: 



^k,k' — -Dk,k' (2.2) 



where 



I^k,k = ^ E (l^(k/2 - q)|'ri,,t,q + l^(k/2 + q)| Vi,q) • (2-3) 

O q 



Here Uh^a is the number operator for holes (1 — nc,a where nc,a is the number operator for 
electrons). Since the magnitude of -Dk,k is proportional to density of holes in the antifer- 
romagnet, -Dk,k is smaller than x, the total number of holes in the system. Moreover, the 
largest fraction of holes is actually residing in the stripes. Since ( p.3|) is much smaller than 
one, the violation of the bosonic commutation relations can be disregarded and we can treat 
the bound states as real bosons. At this point it is convenient to count the number, Ne, of 
fermions in the system. Assuming that there are Ns/N stripes of size A^^ (in lattice units) 
the number of antiferromagnetic sites is (A^^ ~ Ng/N) x A^^ and the total number of fermions 
is: 

AT, = N^{1 - 1/N) -2Y^P\k)P{k) + E4,>k,. (2.4) 

k k,(7 

where the factor of 2 comes from the composite nature of the bosons and the last term is 
the number of fermions in the stripes (see below). 

The Hamiltonian of the composite bosons can be written as 

Hp = E(-Ek + 2^)ptp^ + ^ E (2-5) 
k ^ i 

where i?k is the dispersion relation of the bosons and U is the local repulsion between pairs. 
In ( p.5|) fi is the chemical potential for the electrons and it comes with a negative sign 
because the chemical energy appears as —fiN with the total number of electrons given in 

(ED- 

Suppose there are no other interactions in the problem (that is, stripe fermions and d^2_y2 
bosons do not interact). The Hamiltonian (|2.5| ) is the so-called Bose-Hubbard model. This 
problem has been studied to a great level of detail and its phase diagram is well-knownil. 
Let us consider the situation close to the minimum of i?k at, say, k = K. Expanding close 
to this point one can rewrite the Hamiltonian for the bosons as 



H = Y. (-Ek + 2fi)N, + ^N^] -IbY. P}P^ + h.c. (2.6) 



where Ib is the hopping energy of the boson. The phase diagram is strongly dependent on 
the boson chemical potential, /^b, that is given by 
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(2.7) 



In the absence of disorder the superfluid phases are separated from the Mott insulator 
phases by hnes of second order phase transition. At = the Mott insulator phases 
extend in the range n — 1 < fis/U < n where n is a positive integer that gives the number 
of bosons per site. Notice that because of condition ( |1.20|) we must always have /xs < 0. 
Therefore, the only allowed state for the bosons is the Mott insulator with n = 0, that is, 
the vacuum. For a finite value of there is a critical value of the hopping above which the 
system becomes a superfluid. However, even a small amount of disorder suppresses the Mott 
insulator-superfiuid transition by the creation of an insulating state. Disorder is unavoidable 
in these systems since the charge comes from counter-ions out of the Cu02 planes. Thus, 
the conclusion is that if the bosons are decoupled from the stripes the Bose system is in the 
insulating Bose glass state. 



III. STRIPE FERMIONS 

As we have argued the single particle excitations at the regions of high density of the 
CCS is due to the magnetic confinementS. In the absence of the CCS the holes can move 
freely with hopping energy t in both in the x and y directions and the Brillouin zone is 
defined for —n/a < kx,ky < it /a. In the presence of the CCS the translation symmetry 
of the lattice perpendicular to the ordering vector is broken. Suppose that the stripes are 
separated by a distance Na. Then the new Brillouin zone is given by — 7r/a < < n/a and 
—T[/{Na) < ky < Ti/{Na). Thus the original band has to be folded back into this zone and 
gaps open in the single particle spectrum generating new bands. 

When N ^ 1 the tunneling between stripes is suppressed because of the large distance. 
A WKB estimate of the transverse tunneling energy using the confining potential ( |1.9| ) gives: 

^ ^"^^!rr^' e-^^^^"V^^^ (3.1) 



where m is the hole mass. Observe that for A^ ^ 1 (|1.18D shows that ps ~ and 



therefore the stripe fermions have large gaps in their spectrum. Moreover, is exponentially 
suppressed at large A^ and intra-stripe charge fluctuations are suppressed. As A^ is decreased 
the confining potential gets weaker and eventually vanishes at A^^ given in ( |1.19|) . Close 



to Nm the system is essentially 2D and the confining potential Vciy) can be treated as a 
perturbation of the atomic one and small gaps of order ps{N) will open at 7r/{Na). At 
this point, A^ N~, the bandwidth of the lowest band is simply 4t(l — cos(7r/A'rn.)) < 4t. 
Exactly at A^ the stiffness vanishes and the system becomes 2D since the holes can move 
freely, that is t_L(A^ > A^) = t where t is the hopping energy in the absence of the CCS. 
A simple and convenient way to parameterize the hopping in the transverse direction in the 
whole parameter range is 

t^{N) = (3.2) 

where Aq ~ A"™/ ln[l/(l — cos(7r/Am))]. The hoppine along the stripes, is not so sensitive 
to the distance between stripes. It has been shownli3 that for A^ — oo the hopping energy 
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is reduced from t (as in the uniform system) due to the dressing by strings to a value of 
the order of J ~ t. On the other hand when 1 one has t|| = t. Thus, t|| is a smooth 
function and its variation with can be ignored. The anisotropy in the hopping energies 
will lead to anisotropics in the Fermi surface as well. Suppose the stripes are oriented along 
the X axis. The stripe fermion dispersion relation can be written as 

ek(A^) = -2t|| cos(A;^.a) - 2t^{N) cos{Nkya) (3.3) 

and it is strongly anisotropic as long as tn ^ t± (this is always true for > Nm)- In the limit 
of ^ 1 the transverse component can be completely disregarded and the Fermi surface is 
obtained by filling up all the states up to kx = ^kp where kp is the Fermi momentum (see 
Fig.^(a)). Moreover, the Fermi surface is open along ky since the system has no dispersion 
in that direction. If the stripes are oriented along the diagonals on a square lattice then in 
momentum space the dispersion looks like the one given in Fig.^(b). Then it is convenient 
to rotate the orientation of the lattice by 7r/4 and to change the lattice spacing from a to 
a/\/2. Introduction of t± leads to a change in the curvature of the Fermi surface by an 
amount 6k± close to the Fermi points. 

As it has been noted in the context of organic superconductorsH the strong anisotropy 
of the hopping produces strong dependence of the physical properties with temperature. 
At high temperatures, that is, T > t_i_{N), the uncertainty on the transverse momentum is 
larger than the size of the Brillouin zone, Sk_i_ > 7r/{Na), and the curvature of the Fermi 
surface is blurred by thermal effects. In this case the motion of the electrons perpendicular 
to the stripes becomes incoherent and electrons become confined by thermal effects to a 
region of size Na. Thus, the motion of the electrons is essentially ID and the 2D aspect of 
the stripe problem can be disregarded. In the opposite limit, T < t_L(A^), the curvature of 
the Fermi surface is larger than the thermal effects and 2D physics start to play a role. So 
there is a ID to 2D crossover as a function of temperature that occurs at Tcr{N) ^ tj_{N). 

In most of the paper we focus in the region where A^ < A^^ so it makes sense to talk about 
stripes and the existence of a CCS. In this region the transverse hopping given by ( ^.2[ ) is 
exponentially small and the problem can be treated as purely ID. Consider a single isolated 
static a stripe that we describe as a ID electron gas with dispersion e^. The non-interacting 
Hamiltonian is simply 

Ho = Y.i^k-f^)^pl^^Pk,a. (3.4) 

k,a 

As usual in the case of the ID system we will assume that the physics is dominated by the 
excitations close to the Fermi points. Expanding the dispersion close to ±kp {e±kp = we 
can write the above Hamiltonian as 

Hq = ivFj2 / {'^RAs)dsipRAs) - ^I,a(s)5s^L,a(s)) (3.5) 

where s is the distance along the stripe (if the stripe is oriented along the x axis then s = x, 
but if the stripe is diagonal then s = a/x^ + y^) vp is the Fermi velocity and 'ipnis) and 
'4^l{s) are of right and left moving electrons, respectively. The original fermion operator is 
written as 
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Ms) = ^rAsY''^' + ^L,<x(s)e-^^^^ (3.6) 
and these can be bosonized via the transformationii 

The bosonic modes can be described in terms of amphtude, and phase, 6^, bosonic 
modes as (Pr^lA^) = 0o-(s) =F Gu{.s). The bosonic fields can be rewritten in terms of charge 
and spin bosonic modes: 0c,s = (0t ^ 'Pi) and 6c,s — ^{^i ^ ^i)- the average density 
on the stripe is po then we can write, Scr=t,j, '^i('5)V'o-('5) = a(po + ^p) and 5p = ds4>p/ y/n in 
bosonized form. The Hamiltonian of the problem for forward scattering readJ^l: 

Hf= E^-^ [ds {k, {d,e,f + K~^' {d.ct^.f} (3.7) 

where are the Luttinger parameters (dependent on the electron-electron interactions) 
and 

= ^ (3.8) 

are the charge and spin velocity. For the non-interacting system we have Kp = Kg = \. By 
the same token the action that describes the bosonized fermions is given by, 



S = T.7r^ fds fdr [{dre,)' + vl [8,6,)' 

= E jdsjdr [{dA,f + vl [dgc^.f] . (3.9) 

Besides forward scattering we may also include a backscattering term: 

Hi(x j ds ^lj,t(s)V'ia(^)^-R.i(^)^^.T(^) + h.c. 

= gi f ds cos(v^0s(s)) (3.10) 



where gi is the coupling constant. In principle a Luttinger liquid can also include Umklapp 
terms. The Umklapp terms are only important at commensurate fillings and especially at 
half-filling when they are responsible for the Mott-Hubbard gap. However, half-filling for the 
stripes implies that there is one electron per site, that is, the stripe is depopulated. This can 
only happen if the boson level crosses the chemical potential and the stripe state becomes 
unstable. Throughout the paper we assume that the stripe state is stable and therefore away 
from half-filling. Experiments by Tranquada and collaborators in LNSCO indicate that the 
stripe filling is 1/4^3. Thus, we disregard the Umklapp processes. 

The importance of the backscattering term can be readily understood by a simple per- 
turbative renormalization group (RG) calculation. It is a simple exercise to show that if we 
shrink the bandwidth of the stripe electrons from Ec to Ec — dE^ the coupling constant gi 
renormalizes according toii 
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dgi 
di 



2{1 -Ks)gi{i) 



(3.11) 



where i = ln{l/ Ec). Thus, gi is irrelevant for < 1 and relevant for Ks > 1. For repulsive 
interactions it turns out that Ks < 1 and therefore the interaction is irrelevant and the 
system is described by a gapless Luttinger hquid. If interactions are attractive Ks > I and 
the backscattering is relevant leading to the opening of a spin gap in the system signaling 
tendency to superconducting fluctuations. We have started from a repulsive model and 
therefore we assume all the interactions to be repulsive throughout the paper. Although 
the backscattering term is irrelevant it is well-known that the slowest decaying correlation 
function for a repulsive interacting Luttinger liquid is the CDW one (it is equivalent to the 
case of Kp < 1). 

If the CDW is commensurate with the lattice (as it seems to be in the case of LSCO) then 
the lattice pinning is strong and the system becomes insulating^. If the CDW is incommen- 
surate then a Luttinger liquid state is possible but any amount of disorder will lead to CDW 
pinning driving the system again towards an insulating regime. This is certainly undesirable 
from the phenomenological point of view since these materials are still metallic at finite 
temperatures. Many approaches that start from a purely ID description of the stripes are 
confronted with this serious problemE^I. These approaches usually invoke the introduction of 
long-range interactions among the stripes. It is hard to reconcile the long range interactions 
between different stripes with the short range interacting within the stripes. For that to 
happen the dielectric function of the material should be extremelly anisotropic ( experimen- 
tally it is known that this is not the cas We believe that once the long-range forces have 
done their work in creating the stripe array the final metallic state screens any long-range 
forces left, leading the system to purely local interactions. If the inter-stripe interactions 
remain long-range the ground state should be a Wigner crystal^ instead of a CDW and in 
both cases the stripes should be insulating. We stick with the assumption that after the 
CCS is formed the interactions between the elementary excitations discussed here are local. 
The ground state is a CDW state and therefore insulating. In fact this is the picture that 
arises from mean-field approaches to the stripe problem since in these mean-field studies the 
stripes are non-interactingil'ii. 

In the case of the CCS state where there are an infinite number of stripes we add a new 
index to the fermion operator, ipn,a{s), where n labels the particular stripe. The Fourier 
transform of the electron operator can be defined as 



where the anti-commutation relations are preserved {{ipn,a{s),ipm,a'{x')} = S{x — x')6n,mSa,a' 



and {^.,k,V'a',k'} = (27r)25(k-k')). 

IV. COUPLING BETWEEN STRIPE FERMIONS AND 2D BOSONS 



As we have discussed in the last two sections in the absence of coupling between stripe 
fermions and bosons the ground state is insulating: a Bose-Mott insulator in the antiferro- 
magnetic regions and a CDW state along the stripes. In order to determine the coupling 



nir/a fj^j^ n 

J-ir/a ZTT J- 



/-/(^«) dky 

-TT/{Na) 2tT 



(3.12) 
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between bosons and fermions we consider simple arguments based on conservation of mo- 
mentum and energy. 

As we have argued, the nature of the elementary excitations in the problem depends 
strongly on their position in real space. On the one hand, bound pair bosons are not 
eigenstates of the Luttinger liquid. On the other hand, Luttinger liquid bosons (collective 
sound waves) are not elementary excitations of the doped Mott insulator. Thus, the simplest 
process that couples these two types of excitations is a decay process in which bosons are 
continuously transformed in stripe fermions and vice- versa. 

Let us consider the case of static stripes. The relevant coupling that conserves linear 
momentum can be written as: 

^k(q)P-k^q,TV'k-qa , (4-1) 

where Vk(q) is the boson-stripe coupling that, by conservation of angular momentum, has 
dx2_y2 symmetry (that is, Vk(0) oc cos{kya) — cos^kxo)). This is a four-fermion coupling 
in terms of the original operators (see ( |2.1| )) and describes a fermionic scattering process 
where in the final state two electrons always end up in a bound state. This kind of coupling 
looks unusual because there are only annihilation operators and it corresponds to the vertex 
shown in Fig.|](a). It only happens this way because we use the electrons as assymptotic 
states in the scattering process. In fact, (|4.1| ) is similar to fermion-boson models that have 
appeared in the literature of superconductivity over the yearsHS Notice that since the 
boson is composed of holes the direction of the current is inverted (holes moving forward 
are equivalent to electrons moving backward and vice- versa). Thus to destroy a boson is 
equivalent to create two electrons in the antiferromagnetic background. In the process ( [4.1|) 
a pair of holes is destroyed in the antiferromagnet while a pair of electrons from the stripe 
hop into the antiferromagnetic media (see Fig.|^(b)). 

Let us assume that the coupling between bosons and stripes is weak so that the important 
degrees of freedom only exist close to the chemical potential. Thus, when a boson hits a 
stripe there are two main process for decaying. In the first one the two stripe fermions have 
opposite momentum, that is, q ~ kp^ and k — q ^ —kp^ (as shown in Fig.^(b)). Thus, 
linear momentum conservation implies k = 0. Moreover, angular momentum cannot be 
conserved in the scattering process because, relative to the stripe, the two outgoing holes 
have zero angular momentum while the pair have finite angular momentum. Therefore, we 
need a sink of momentum. The other process corresponds to boson decay into two fermion 
moving in the same direction, that is, q ~ kp:K and k — q ^ kp^K. This process implies 
kx = 2kp and fc^ = 0. This is equivalent to a boson moving in the antiferromagnetic ladder 
parallel to the stripe. Since the boson wave-function is restricted to the antiferromagnetic 
region the amplitude of this coupling is suppressed. The conclusion from this argument is 
that for a static stripe the decay processes have vanishing phase space. 

Consider also the difference between longitudinal and diagonal stripes as shown in Fig.|^. 
For a longitudinal stripe the coupling between bosons and stripes is possible because the 
stripes are oriented along the lobes of the boson dx'^-y'^ wavefunction while in the diagonal 
case the stripes are oriented along the nodes. Thus, the coupling between diagonal stripes is 
further reduced as compared to the longitudinal case. Therefore, for diagonal stripes super- 
conductivity is not possible in our model because of the nature of the boson wave-function. 
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In fact, recent neutron scattering experiments observe diagonal stripes in the insulating 
phase of LSCO while vertical and horizontal stripes are observed in the superconducting 
phase0. Moreover, notice that for a given orientation of the boson wavefunction, horizontal 
(vertical) stripes couple to the positive (negative) lobe of the wavefunction. 

As we pointed out above we need a sink (a particle or collective mode) for the momentum 
transferred from the bosons. Moreover, since the bosons are spinless the particle that carries 
the momentum cannot be a spin-wave or paramagnon. The simplest choice is a vibrational 
mode: the boson collides with the stripe, breaks into two fermions and produces a vibration. 
This vibration is due to the fluctuations of the smectic phase of the stripe array. Again, as 
we stressed previously, stripes are not isolated objects that are independent of each other. 
Long-range forces, like the Coulomb, the Casimir or simply entropic repulsion keeps them 
apart and generate a finite stiffness. 

Consider the coupling of the bosons with the stripe fermions when the stripe fluctuates 
a distance u from its equilibrium position. We write the local coupling as 

P{r)ilj^{r)ilJi{r)6{y - nNa - u{x, nNa)) (4.2) 

where we ensure that the coupling occurs at the position of the stripes. If we assume that 
CCS fluctuations are small compared to the periodicity in the system, that is, \u\ <^ Na, 
the stripes fluctuate in a distance scale much smaller than the inter-stripe distance. We 
can expand the Dirac delta function. The first term is just the static stripe problem we 
discussed above and, as we argued, has no phase space. We therefore keep the next term in 
the expansion leading to a coupling of the form 

- Vk(q)Ap,„P„k^q,T^k-p-q,i (fep,a + &-p,a) (4-3) 



where Ap^ (x Py ^ Q when p ^ is the stripe- vibration coupling constant. The two 
processes associated with the interaction in ( |4.3| ) are shown in Fig.|^. 

Before we study ( [4. 3] ) in detail it is worth investigating the problem via second order time 
dependent perturbation theory. The transition probability between non-perturbed states in 
the system is given by: 

W{t)^ E A\VM)?\\A\l - n.){l - Tr.^p-q) '"7'^V^i^"° ^^^^^ (4 

where rlq is the Fermi-Dirac distribution function. Observe that the probability amplitude 
only grows with time for states such that 

- -^k ± ^p,a = -^qx - ^kx+Px-qx (4-5) 

where the plus (minus) sign implies absorption (emission) of a vibration. Remembering 
that eg ~ /i we obtain the condition that —Eo^ky ± ^o,ky,a = — 2/i. Moreover, because 
vanishes at /c^ — »• the largest coupling between bosons and stripes happens at finite 
momentum. Indeed, the coupling increases as one goes towards the Brillouin zone edge. 
Thus, the largest coupling in the problem occurs at K = (0,7r/a), that is, perpendicular to 
the stripes. It may be somewhat surprising that the most important part of the physics 
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occurs in a direction perpendicular to the orientation of the stripes but this effect has been 
already observed in numerical simulations of spin-fermion models with stripe formation!. 
This result indicates that in the case of the boson condensation we should expect it to occur 
close at K (as photo-emission experiments indicateil) but not only that: we expect a CCS 
distortion at the same wave-vector so that the momentum carried by the superfiuid bosons 
is transfered to the lattice. Thus, there can be a double condensation: the superfiuid density 
is modulated with finite wave- vector K and the CCS deforms in order to follow the variation 
of the superfiuid density. Therefore a condensation of the CCS at finite wave-vector will 
induce phonon anomalies at the same wave-vector. Thus, we expect phonon softening close 
to (0, 7r/a) together with the condensation of bosons in the same point in momentum space. 
In the normal state the situation is depicted in Fig.||: there is no superfluidity and no lattice 
distortions. In the ordered phase we expect the situation shown in Fig.^ where the superfiuid 
density is modulated with wavelength a/vr in the direction perpendicular to the stripes and 
the lattice is distorted in the same direction with the same wave-vector. Notice that there is 
a doubling of the unit cell in the direction perpendicular to the stripes and therefore the CCS 
is dimerized. This effect is similar to the Peierls distortion in ID systems or the Jahn- Teller 
effect where there is a lowering of the symmetry with a simultaneous gain in energy!!. In 
momentum space the situation is depicted in Fig.^: phonon anomalies and superfluidity 
should appear at same point on the edge of the Brillouin zone. 

Observe that due to the conservation of energy, — -Eo,fcj; +2yU = ±f2o,fcy,a, and the stability 
condition ( |1.20|) is satisfied only if emission is allowed. One can generate a boson by creating 
two electrons at the stripe at the cost of absorption of a vibration. This can be clearly seen 
in the diagram in Fig.|l| where in order to create a boson two electrons occupying the boson 
state have to be excited to the Fermi surface by the absorption of a vibration. Conversely, 
two holes from the stripe hop into the antiferromagnetic ladder and form a boson. Thus, 
our conclusion is that the most relevant coupling in this problem is given by 

Hc= - Vk(q)Ap,„Pk^q,t?/'_k+p-q,i^{,,a + ^.C. . (4.6) 

k,p,q,a 



V. MEAN-FIELD THEORY 

In order to understand the phase diagram of this problem it is illuminating to consider 
the mean-field solution. Let us rewrite the full Hamiltonian by collecting all the terms given 
in the previous sections (we neglect the interactions between electrons in the stripe for the 
moment being). The Hamiltonian is written as 

H = Y: [i-E^ + 2/x)P^Pk + n^bib^] + ^(ek - /x)<,^k,. 

k k,(7 

+ I E - E ^k(q)ApPk^q,TV-k+p-qa&i. + h.c. (5.1) 

i k,p,q 

where we have simplified the problem by considering a single polarization of the vibration 
field perpendicular to the stripes. 
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Notice that the Hamiltonian is invariant under the gauge transformation 



P ^ P e""^^ 
b^b e+*'^^ 

implying the phase of the bosons and vibrations are coupled through the fermionic degrees 
of freedom. This symmetry is the one that is broken in the ordered state. 

In the mean-field theory we split the coupling term into three different pieces leading to 
a mean-field Hamiltonian of the form: Hmf = Hp + Hb + Hs, where: 

Hp = Y: [i-E^ + 2^)PlP^ - gp,^P^ - g*p^^Pl\ + y E 



k 



Hs = I](ek - /i)^k,>k,<7 - J2 9s{K P, q)^/'q,TV'-k+p-q,i + h.c. (5.3) 

k,(T k,p,q 



where 



9pm = E Vk(q)Ap(z^q,tZ^„k+p-q,i) {bp) 
p,q 

9B,k = AkEVp(q)(Pp) (V^q,t^p+k-q,i) 

(75(k,p,q) = rk(q)Ap(Pk)(M (5.4) 



are the mean-field coupling constants. From ( |5.3| ) it becomes obvious the symmetry breaking 
processes in the ordered phase: spontaneous gauge symmetry breaking leading to superflu- 
idity of the bosons, superconductivity of the fermions and spontaneous symmetry breaking 
of translational symmetry of the CCS. 

Thus, for a transition to a non-trivial conducting state the coupling AV^ is fundamental. 
The simplest conducting state is the superfiuid state where the charged bosons condense into 
a macroscopic condensate. In accordance with the discussion in the previous section, conser- 
vation rules tell us that the largest coupling occurs at the largest momentum perpendicular 
to the stripes. Suppose the bosons condense at k = K, that is, (Pk) = {Pk) = V^b- In 
order to conserve momentum in the interaction in ( p.l] ) and preserve time reversal symmetry 
one should also require the vibrations to condense, that is, (bj^) = (6k) = y/Ny. Notice 
that the Cooper pairs have zero center of mass momentum. This is equivalent to a static 
distortion of the lattice with modulation K. Thus, at the mean-field level, a transition to 
a superfiuid state is accompanied by a static lattice distortion. We notice here that the 
choice of the largest coupling as the one that drives the system into the ordered state is the 
traditional one in mean-field theories. In fact, the doubling the unit cell proposed here also 
occurs in some large N treatments of the t — J models without phononsHS. Thus, this 
tendency of lowering the energy of the system by breaking the lattice symmetry seems to 
be quite generic. 

The solution of the bosonic problem is straightforward because it only requires a shift of 
the boson operators by a constant: 
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(^K 



-Ek + 2fx 

gB,K 
n 



(5.5) 



K 



implying that the actual CCS distortion is given by 

V2^gB,K] 



6u{K) 



(5.6) 



where My is the mass associated with the vibrations. 

The electronic problem can be also readily solved since it is a problem of electrons with 
a pairing potential. We can write Hs as 



hs=j: 



ck - ^) I] V^k,>k,a - AkV^k,TV^-ka + h.c. 



where the superconducting gap is given by 



Ak = AKaV^Vy ^k(O) 



(5.7) 



(5.8) 



and therefore has dx^-y^ symmetry as the pairs that generate it in the first place. Here we 
have defined 
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V 



S 



as the superfluid and distortion densities, respectively. The Hamiltonian 
onalized by a Bogoliubov transformation: 

V'k.T = cos(6'k)rfk,T + sin(6'k)4a 
^-k,T = cos(6'k)(ik,i - sin(6'k)c?k,T 



where 



tan(2^i, 



ek - /i 



and 



Hs = Yl ek4,kC^<T,k 



(5.9) 
can be diag- 

(5.10) 
(5.11) 
(5.12) 



where 



Ck = V(^k - /i)^ + |Ak| 



(5.13) 



is the quasiparticle dispersion. 
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The free energy per unit of area, f{T), in the system can be obtained in a straightforward 
way: 

/(T) = (-Ek + 2/i)a° + )2 + Qj^a'y + W [(e^ - fi) - 6^(1 - 2n^)] (5.14) 



where 



2 S ^ 



1 

= , , (5.15) 



e/3ek + 1 



is the quasiparticle occupation. We minimize ( 5.14 ) keeping the number of particles fixed. 
We obtain two equations: 

2b ^ Ck 

= E ^tanh{Pe^/2) (5.16) 



and therefore 



n^a^y = al{-E^ + 2^i + Ual) 

^U{air (5.17) 



that is, the superfluid and vibration density are tied to each other. The last line in ( ^.17| ) 



comes from the fact that we expect U ^ — -Ek + 2/i. Using the above equation we reduce 
the mean-field problem to a single equation 

\2 2 ta2 

^^K = E —tanhiPe^/2) . (5.18) 

2o Ck 

Notice that although the coupling constant has d-wave symmetry the Fermi surface of 
an infinite array of ID stripes has no curvature (as shown in Fig.^(a)) because the fermions 
do not propagate in the direction perpendicular to the stripes. Thus, although the order 
parameter in (|5.8|) has d^.2_y2 symmetry, the Fermi surface is fully gapped. The bosons have 
a tendency to produce a d-wave superconductor but the condensate only takes advantage of 
one of the lobes of the boson wave-function. Thus, the final result is an s-wave state. As we 
will discuss in section 11^ a true d-wave superconductor can only develop if the system is made 



out of domains with horizontal and vertical stripes or a one particle tunneling is included 
(changing the fermion dispersion relation as in ( |3.3| )) and the Fermi surface becomes rounded 
leaving space for a true d-wave state (although rather anisotropic). This seems to indeed be 
the case of these systems where microscopic twin boundaries are unavoidable. Because the 
domains are large the superconducting properties of the domains can be studied as if the 
domains are macroscopic superconductors. Taking the geometry of the Fermi surface into 
account we can simplify the above equation by using the fact that ( |5.15| ) only depends on 
kx and therefore the sum over ky just gives Ng/N leading to a simpler equation 
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JM Jo /TT Cfc 

where Vq is the strength of the couphng Vk(q) at the stripe and = 



Cfc — /u)^ + A2(cr^)2(— i?K + 2/i + ?7o"^)/^^K- This equation can be solved numerically 
assuming e^^ = — t cos(A;a;a). 

We can now study the phase diagram as a function of the coupling XV and the dis- 
tance between stripes A^. We expect the characteristic energy scales in the problem to be 
dependent on the distance between stripes. This is certainly true for the case of antifer- 
romagnetic order where the distance between the stripes determines the strength of the 
quantum fluctuationJll. Thus we expect the phonon and boson spectra to change with A^, 
that is, there must be a continous change in the spectral weight with the interstripe dis- 
tance. In this paper, for simplicity, we will assume that this change is small and the density 
of states does not change with N (see, however,!). In Fig.|Tl| we plot function of 



for Ek = 1.9t, /i = -0.6t, Qj^ = 1.2 x IQ-H, U = it, \V = 0.5t at T = 0. Notice 
that there is a discontinuous jump in cr^ at N^p ^ 6 and it grows in a quasi-linear way 
with 1/A^. In Fig.[l^ we plot cr^ as a function of T/t for the same parameters given above. 
Notice that cr^ is weakly dependent on the temperature up to T = Tgp ~ 3.5t where it has a 
discontinuous jump to zero. As one can see from the numerical solution Tgp is of the order of 
the bandwidth of the stripe fermions and therefore very large. As we show below this is not 
the actual transition temperature because we have not included phase fluctuations into the 
problem. Thus Tgp (like T* defined in Section |) is a crossover temperature scale. We note 
that Tgp ^ T*, that is the mean-field critical temperature is larger than the temperature 
scale above which makes sense to talk about boson pairs and confinement. Therefore, Tsp is 
actually not observable since the theory brakes down at T*. 

In order to understand the nature of the phase diagram we simplify the problem by the 
introduction of a cut-off A of the order of the inverse of the lattice spacing and linearize the 
spectrum at the chemical potential. The integral simplifies and at T = reduces to: 

271VfN I Ao ' 



where Aq = VoXa%y {—Ej^ + 2yU + Ua'^)/VlK gives the amplitude of the gap. Notice that 
the transition at the mean-field level is of first order because the r.h.s. of ( |5.2CI| ) vanishes 



as CTjj in contrast with the usual BCS equation that contains only the logarithm and 
is finite in this limit. Some analytical progress can be made in the limit of large U. The 
critical value of A^ and the value of the superfiuid fraction are given by: 



A^ 



ln(l + V2) (AV-o)^/3 

2 (7rt;^/a)i/3l]^/3t/i/3 



that determines the zero temperature mean-field phase diagram. As we vary A^ there is a 



critical value of Vq (that we call Vc) below which order is not possible. In Fig.O we plot Vc 
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as a function of We find that Vc oc A^" with a > 0. This result demonstrates that as 

the stripes get further apart a larger coupling is required to stabilize long-range order in the 
system. 

As pointed out by EK, the fact that the superfluid density is so low implies large phase 
fluctuations^. If one includes the phase of the order parameters (as in (|5.2| )) long-range 
order requires that 

(PK) = (|PK|)(e-^^-)^0 

(&K) = (|6K|)(e+^^n^0 (5.22) 

where (|-Pk|) oc \fo^ and (|&k|) oc \fo^ are the amplitudes of the order parameters. Thus, 
ordering also requires (e^*"^") 7^ with a = B,V. In our mean-field theory we have fixed 
arbitrarily ip^ = 0. This is not correct since fluctuations of (pa, especially in the 2D system, 
are fundamental for the development of long-range order. In the next sections we discuss the 
role of these fluctuations in determining the actual phase diagram. In other words, at high 
temperatures bosons and superconducting stripe fermions are formed but they are incoherent 
because of the existence of topological defects: vortex-antivortex pairs in the superfluid and 
dislocation loops of the distorted CCS. The interesting result of our calculation is that 
although the amplitudes of the two order parameters appear at the same in temperature T* 
(as given by (|5.17|) ) the actual ordering temperatures can be different because they depend 
on the phase stiffness of the relevant degrees of freedom. 

VI. PHASE FLUCTUATIONS 

As we discussed at the end of the previous section the problem of phase fluctuations is 
fundamental for the description of the ordered state. Therefore the topological excitations 
will be fundamental for the determination of the properties of the superconducting statei^. 
The simplest way to discuss the importance of phase fluctuations is by studying the partition 
function of the problem that can be written as: 

Z = J DPDPD^D^DlDh e- ^0 dTL[p,p,^,i.rbM 

where 

L = j dh \p{r){dr + ^(V) + 2/i)P(r) + ^N^{t) + h{Y){dr + w(V))6(r) 

+ J dx J dx''^llJn^„{x)tjjn,a{x)W^y{x - X')llj^y{x')'lljn,a'{x') 
a, a' 

~E/ dxV{x^nNa) ^P{x,nNa)'ipn,^{x)'ipn,i{x)h{x^nNa) + c.c)^ (6.2) 



is the Lagrangian associated with Hamiltonian ( |5.1| ) and j3 is the inverse temperature. Here 
we have redefined AV^ V and introduced the electron-electron interactions W^^^i {x — x') . 
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We redefine the fields in order to separate the problem into slow and fast variables: 

1 

71 



P(r) 
6(r) 



1 



V2 



h+{r)e 



iK r 



iKr 



(6.3) 



and perform the same expansion as in the case of fermions (see ( |3.(j| )). 

Substituting (|3l6|) and ( |6.3|) into the Lagrangian ( |6.2|) we find many terms that oscillate 
fast. We disregards all the oscillating terms and study slow modes only. The Lagrangian 
can be written as 



a=± n 



(6.4) 



where 



2Mf 



Pair) + ^N\v) 



UJQ 



2M, 



&a(r) 



(6.5) 



V . 



where we have assumed that K is a local minimum in the energv of the bosons and vibrations. 
This assumption is supported by studies of the t — J modeled and from the fact that the 
vibration mode is at the edge of the Brillouin zone. Here, Mb is the effective boson mass. 
The form of the non-interacting boson-vibration Lagrangian given above is required by 
Galilean invariance of the problem in the long wavelength limit. Although anisotropics in 
the dispersion of the modes close to K should exist, we disregard them since they are not 
fundamental in our discussion. Furthermore, 



(x) {dr + iVpO:,) i'n,RA^) + '^'u^lA^) (^^ " '^'"pdx) 1pn,LA^) 



+ J2 [Wk=o{^n,RAA^^,RAA + ^n,LAA'^ri,L,a{x))x 

X ii>R,a'ix^n)^pRy{x,n) +^i_^,(x,n)^L,<^/(x,n)) 
X {^R,a'ix,n)^L,a'{x,n) + tljL^^,{x,n)iljRy{x,n))^} 



(6.6) 



is the isolated stripe Lagrangian [W^ is the Fourier transform of the electron-electron inter- 
action) and 



= -^Yl J dxV{x,n) P+{x,n)b+{x,n) + P^{x,n)b-{x,n) [ipn,R,^{x)ilJn,LAA 



+ '4^n,L,l{Ai'n,R,'\{x)]+ C.C.} 

is the boson-stripe-vibration coupling. 



(6.7) 
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Due to the symmetry we simplify the problem by changing to amplitude-phase modes 



via: 



P^ = P- 



-tipB 



and use the bosonization technique described previously so that 



(6.8) 



(6.9) 



a=By 



where LiD[On,p, 4>n,s] is given in ( ^lOf ). Furthermore, 



(6.10) 



with Eb = EiQ, Ey = f^K, Ub = U and f/y = and, 

dx , 



Lc = - V / — \/(a;,niVa)y^i^e-'(^«("'"^'')+^^("'"^'^)-^'-''(")) cos(v^0„,,(a;)) + c.c. . (6 
„ J Tia 

Observe that the gauge symmetry discussed in the previous section is explicit in the 
Lagrangean. In order to gauge away the phase fields we define a new bosonic field 9c- 



On,c{x) = 9n,p{x) - 



^b{s^i nNa) + '{'vi.x, nNa) 



2ti 



(6.12) 



so that Lc simplifies to 

dx 



/djr 
— V{x, nNa)y^o^ay cos{V2nen,c{x)) cos{V2n(j)n,s{x)) (6.13) 

and in the final Lagrangian we replace 6p by 6c and add a new term that reads: 



(<9r(v5B(x, nNa) + ipv{x, nNa)))' 



+ vl {dr,{ipB{x, nNa) + (pv{x, nNa))y 



+ 



27TV, 



[dri^six, nNa) + (fvix, nNa))dr9n,c{x) 



+ vldx{LpB{x,nNa) + ipv{x,nNa))d:cOn,c{^) } 



(6.14) 



The physical interpretation of (|6.12|) is quite interesting. Observe that the fermion current 
along the stripe is 



jn{x) = —=d^en,p{x) 



-^d^en,c{x) + -^d^LpB{x,nNa) + -^d^Lpv{x,nNa) 
yT^ 7rv2 7rv2 



(6.15) 
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where the first term is just the normal current, the second term is the superfluid Josephson 
current and the last term is a current that is driven by the dislocations in the CCS. Moreover, 
observe that the bosonic field 9c is still coupled to the phase fields through ( |6.14| ). Thus, 
the stripes carry along both superfiuid and CCS distortions. 

In order to proceed with the calculations we study the phase fiuctuations around the 
saddle point equations for the complete action. At the saddle point we have: cTq, = cr°, 
(fa = with a = B,V. This just gives the mean-field solution given in Section |V] with the 
difference that at the saddle point the stripe action reads (we drop the subscript n): 



S 



/3 



drdx i — - 
2vo 



+ 



{dr(t)s{x, T)f + vl{d^(t)s{x, T)f 



/3 



— j dr I dxV cos{\''2n6c{x)) cos(v 27r0s(a;)) 



(6.16) 



Eq. ( |6.16| ) is the interacting version of the Hamiltonian used in the previous Section with 
V 



/a%ayV. The last term in ( |6.16|) is the Josephson coupling induced between different 
stripes by the bosons. A similar coupling was proposed for the case where stripes cross 
each otheiS. The effect of the pairing term can be easily seen in a renormalization group 
calculation in first order in V. Repeating the RG calculation done in Section 
easily show that the dimensionless coupling constant v = V/Ec renormalizes asE 



we can 



dv 




p J \ 



vii) 



(6.17) 



Notice that the operator associated with V is relevant when Kg + K^^ < 4 and irrelevant 
otherwise. For a Hubbard model with local repulsion Uh it requires that Uh > l-8t for the 
operator to become irrelevant. Thus, even if the stripes are deep inside of a CDW state the 
coupling V drives the system into a superconductor. When V is large then the last term of 
( |6.16| ) acquires a expectation value and the bosonic fields get pinned at their minima: 



9r 



{n + m) 
{n — m) 



(6.18) 



where n and m are integers. Thus, from ( p.l5| ) one sees that the only currents circulating 
along the stripes are the Josephson currents due to the superfiuid and CCS fiuctuations. 
Fluctuations around the minima (|6.18|) are massive (although a superfiuid current still fiows) 
and a spin gap opens in the spectrum. Indeed, a simple integration of ( |6.17D gives 



Vr 

Er 
9 




(6.19) 
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where Vr and Er are the renormahzed couphng constant and bandwidth. Observe that the 
RG flow stops when Vr ~ Er ^ m where 



m 



^'4 



1/(1-9) .y.g/il-g) 



(6.20) 



gives the amphtude of the spin gap (when g = 1 the operator is marginal and m ^ 
Ecexp{—Ec/V}) . Since the stripe modes are massive the coupling terms to the phase 
modes in ( |6.14]) are suppressed. 

The calculation can proceed in the usual way by expanding the action around the fluc- 
tuations of the amplitudes to second order (that is, we write (Jq, = 0"° + 6aa and integrate 
over 6a). The final result reads: 



J2 dr(Pair,T 



+ 



2TiVpNa 



^a=B,V 



Ka=B,V 



r,r 



(6.21) 



where the factor of appears because we have coarse-grained the fields in the direction 
perpendicular to the stripes. 

The parameters that appear in ( |6.21| ) can be obtained directly from the symmetries of 
the original bosonic action ( |6.2| ). The bosons obey periodic boundary conditions {P{r,f3) = 
P(r,0), b{r,j3) = 6(r,0)) while the fermions obey anti-periodic boundary conditions in the 
imaginary time direction {ip{r,P) = —ip{r,0)). Suppose we change the boundary conditions 
so that P(r,/3) = P{r,0) exp{—i6(p}. From ( |6.2|) we can enforce the original boundary 
conditions if we re-interpret the problem as a shift in the chemical potential of the bosons 
from n to n' = fi — Because bosons and fermions are coupled, we are forced to 

impose new boundary conditions for the fermions, namely: ip{r,f3) = —ip{r,0) exp{i6(f/2}. 
The whole action is invariant under the change in the boundary conditions if we assume 
a shift in the chemical potential. If the shift is infinitesimally small we can calculate the 
change in the free energy due to the change in the boundary conditions: 



(6.22) 



and since 5ip/ {3 ^ we immediately find that 



N^a'^ Kb 



— -N 

Q2p 



(6.23) 



where A^e is the average number of electrons in the system. Since the number of holes in the 
system is x we see that = (1 — x)/a?. Thus, we conclude that is the average planar 
density of electrons. 
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Notice that so far this is the first equation where the number of holes appears explicitly. 
In principle the distance between stripes, N, and x should be related but in our theory N is 
an input. When doping is increased the number of holes in each stripe and/or the distance 
between stripes can change. In the normal state the stripes are in a CDW state and there 
is a gap to charge excitations because of commensurability. Since the holes are injected into 
the system at high temperatures during the annealing of the alloy it seems reasonable to 
assume that the CDW state will be charge rigid, that is, instead of increasing the doping of 
individual stripes one would get more stripes but with the same linear doping, rig. If this is 
the case then a simple relation exists between x and N, namely, 

a; = ^ . (6.24) 

The value of Ug is determined by the competition between the gain in the kinetic energy 
along the stripe versus the loss of energy due to the formation of the ADW00. This is 
not a problem we have addressed but DMRG calculationsEl and other approaches estimate 
that the minimization of the energy occurs at = 1/2 as proposed by the experiments 
performed in LNSCOS. Notice that (|1.19|) implies that the pseudo-gap temperature scale, 
T*, vanishes for 1/6 < x = Xgp < 1/4. Assuming that Xgp ~ 0.2 we find that 7 ^ 0.5 and 
there is a reduction of 50% of the antiferromagnetic exchange across the stripe. 
Moreover, is charge compressibility and the superfluid velocity is given by 

CB = . (6.25) 

Kb Mb 

Since the fluctuations of the CCS have zero chemical potential one concludes immediately 
that ay = and is related to the sound velocity in the CCS that is given by 

cv = . (6.26) 

KyMv ^ ' 



It is convenient to rewrite the action ( |6.21[ ) in Fourier space: 
^ = ^Y.\ -o^a^n</'a(k,a;„) + (k^.r^^ + ((1 + A„)A;^ + A;^)) |<^a(k, 

'^P k,n a=By ^ ^'^^ ^ 



K 1 

+ 27rt; iVa '^^"^^^^^^ (<^i?(k, a;^)(^;.(k, ^„) + (^^(k, aj„)v9y(k, ^„))| (6.27) 

where ujn = 2Tm/ (3 is the Matsubara frequency and 

, Kp 

A. = (6.28) 

2-Ka^Na ^ ' 

are the renormalized compressibility and anisotropy introduced by the stripes (we used 
(p.8|)). Notice that both quantities return to unrenormalized values when the distance be- 
tween stripes diverges, that is, iV — > cxd. As a result of the factor Aq, the superconductivity 
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is anisotropic and the correlation lenght, ^, is direction dependent. It is easy to see that the 
ratio of the correlation lengths along the x and y directions is given by 



The key point of ( |6.27|) is that the topological excitations of the superfluid state (the 



vortices) are coupled via the stripe compressibility to the topological excitations of the CCS 
(the dislocations). In a pictorial way consider first the ordered state depicted on Fig.^. A 
defect of the ordered state is shown in Fig.|l4| where a dislocation of the lattice distortion 
leads to a local shift of the superfluid density. In doing so the dislocation can produce 
vortices in the system. The opposite situation is also possible: vortices shift the superfluid 
density and drag the lattice with them. This unusual state of affairs is the result of the 
coupling of the CCS with the lattice and would not happen in ordinary superconductors. 

Although the theory described by ( |6.27| ) is quadratic, the nature of the elementary topo- 
logical excitations is not straightforward and we will leave their discussion for a later pub- 
lication. In this paper we consider a simpler problem of the effective theory for each one of 
the phases. In order to do so we explore the quadratic nature of the action ( p.27| ). 



VII. HIGH TEMPERATURES 

Let us consider the problem at Tc < T < T* so that the amplitude of the order param- 
eters are well developed but phase coherence has not been established. It is obvious from 
( |6.21| ) that this effect is possible due to presence of singular solutions of the field equations, 
that is, due to the presence of vortex-antivortex pairs and CCS dislocations loops. At high 
temperatures we can disregard the time derivatives in ( 6.21 ) since the imaginary time di- 



rection shrinks to zero and the phase fields become independent of r. The effective action 
reads: 

Equation (|7.1| ) describes two 2D XY models coupled along the x direction. 

Let us look first at the effective field theory for v^b by tracing out the ipv modes explicitly. 
This is equivalent to calculate the renormalization of one of the modes by the gaussian 
fluctuations of the other. Thus, in what follows the topological excitations of these fields are 
not directly coupled to each other but only to their "spin- waves" . This approach is valid in 
the weak coupling limit of N ^ 1. It is a simple exercise to show that the effective action 
becomes: 

5B = f E^i^(k)I^B(k)P (7.2) 



where 



n \ 1 2 -^B-Eyk^ + Ec{Eb + Ev)kl 
^^^^^ = ^ Ey\.^ + Eckl ^^-^^ 
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is the effective phase propagator. We have defined 



vf 



c 



2nNa 



(7.4) 



The superfiuid order parameter correlation function can be calculated directly from the 
knowledge of ([7.2[). Indeed, 



(pt(r)P(O)) = as (e 



) 



as exp 



1 



where 



ik r " 



(7.5) 



(7.6) 



is the relevant correlation function. The problem is simplified when we realize that it is 
possible to write 



]_ Ev/Eb 



1 



(7B(k) Eb + Ev k2 Ev + Eb [l + Ec (e^^ + Ey^)] kl + kl 



(7.7) 



that is the sum of the correlation function for an isotropic 2D XY model plus the one for 
an anisotropic 2D XY model. The integrals can be easily done in (|7.6|) and for xA, yK ^ 1 
where A is a ultraviolet cut-off (of the order of the inverse of the lattice spacing) we find 



GB{x,y) 



Ev/E 



B 



In 



+ 



27r/3(Ey + Eb)\II + Ec [E^' + Ey 
1 



A 



\l + Ec {e^' + Ey' 



In 



AJx^ + y 



2jc/3iEv + Eb) 
which, for each direction separately, can be written as 

1 



Gb{s) 

where s can be either x or y and 

Pb = {Eb + Ey^ 



2-K(3pB 



ln(A|s| 



(7.8) 



(7.9) 



Ev/E_ 



B 



1 + Ec {E^' + Ey' 



-1 



(7.10) 



This result can be interpreted as the effective superfiuid stiffness of an isotropic 2D XY 
model. The error in making this approximation is equivalent to an anisotropic change in 
the order parameter (that is irrelevant to the problem of phase coherence discussed here). 
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An analogous calculation can be done for the phase field and it is obvious that one has 
only to exchange Ey by -E^ in the expressions above in order to get: 



Pv = (-E's + E\ 



-1 

EbIE\, 



1 + Ec {e^' + Ey' 



(7.11) 



that can be thought of an effective stiffness for the CCS fluctuations. An interesting conse- 
quence of our calculations is that if My — >• oo then Ey and according to ( |7.11| ) py 0. 



The transition to the static CCS deformed phase is driven to zero temperature. Thus, 
while the system becomes a superfluid-superconductor at T^^b, at any flnite temperature 
dislocations of the CCS do not allow for long-range order. 

We conclude that the transition from the ordered to the disordered phase of the su- 
perfluid or/ and the CCS is due to the unbinding of the topological excitations (Kosterlitz- 
Thouless)liy: vortex-antivortex pairs in the case of the superfluid and dislocation loops of the 
modulated CCS. The transition temperature to the ordered phase can be estimated directly 
from (rrmi ) and (^[TT]) byS 



TKT,a ~ ^Pa (7.12) 



and the superconducting correlation length diverges as 



UT) - aexp , ^ (7.13) 



where 6 is a number of order of unit. 

Although our calculation is completely 2D the transition described here only produces 
quasi-long-range order. True long-range order occurs via the coupling between planes. Thus, 
the real transition into the ordered phase at flnite temperatures is of the 3D-XY type. In 
order to estimate the actual transition temperature, Tc, we assume that the coupling energy 
per unit of lenght, Uc, is small compared to T^t, and therefore the 2D correlation length is 
well-developed when the system undergoes the phase transition. The transition temperture 
is deflned by the amount of energy required to destroy phase coherence between two regions 
of size in different planes separated by a distance c: 

Tc,a^cU,{UTc,a)/af (7.14) 

that is a transcendental equation for for Tc^a- Because of the exponential dependence of ^ 
with the temperature in ( [7.13|) we can solve this equation to logarithmic accuracy: 



This result indicates that Tc depends only weakly on c and it is very close to Tkt- 
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VIII. ZERO TEMPERATURE 



When T — > we can replace the sum over the Matsubara frequencies in ( |6.27[ ) by an 
integral over frequency. Observe that at small frequencies the dominant term in ( |6.27| ) is 
the one that it is linear in the frequency. Thus, as in the case of interactingbosons in 2D 
the universality of the transition is not 2D XY but 2D XY in a magnetic fielcS. We neglect 
all the quadratic terms in the frequency and retain only the linear one. Notice that this is 
equivalent to assume that |ti;| <C o'Blf^B,R < lnVpNaaB / Kp. Thus, the effect of the linear 
term is to introduce a high frequency cut-off given by: 



The equal time correlation function is: 



(8.1) 



(pt(r,0)P(0,0)) = exp^-V ^ 



exp 



-3? 



9b{^) 



E 



1 - e**"'"" - iWBUJ 



9b{\^) 



(8.2) 



where the leading order term in lo vanishes because the integral is symmetric. Observe that 
this the same result of the previous section with the temperature T replaces by 2ujc\ Thus, 
at zero temperature the system will also have unbinding transition of topological excitations 
when N = Nr so that 



TC 



(8.3) 



Usingd]!!), ( lOsD and (|8A|) we see that 



1 - ^ 

2N 



K-B H 



lA) 



2'KVpNa 



is a monotonically decreasing function of On the other hand, from (|7.1CI|) : 



Pb{N) 



(TBiN) 



1 + 



crv{N)MB 
aB(N)/Mv 



M, 



B 



1 + 



ay(N)Mg 
<yB(N)/My 



.5) 



'1+ 



AN) 



that is a rather complicated expression in terms of A^. We have seen from (|5.17| ) that 

CTv{N) 



CTBiN) 



> 1 



:8.6) 



since the electronic energies are much larger than the vibrational ones. Assuming that this 
is the case we can write 
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^^^^^"^i^f + 2vriVaa°(iV) ^^'^^ 

that is a monotonically increasing function of Thus, equation ( |8.3| ) will have a solution 
at l/Ai'c provided that 

Ucil/Nsp) > ^pBil/Nsp) (8.8) 

otherwise there would be no solution and the system would be insulating. Notice that 
Nc < Nsp implying that at zero temperature there there is a region Nc < N < Ngp where 
superfluidity exists even at zero temperature but without phase coherence. Observe that 
according to this theory (even at zero temperature) there must be a region where bosons 
and stripes co-exist but no long-range order is observed. 

Moreover, borrowing the results from ref.S we find that at N ^ Nc we must have the 
superconducting correlation length diverging as 



the superfluid density behaving as 

\[) cx ( 

Nc 



aB{N) oc ( 4 - ) , (8.10) 



and the boson compressibility given by 

where u and z are the critical exponent that at the mean-field level are z = l/u = 2. 
Moreover, close to Nc the transition temperature scales linearly with the superfluid density 

Tc{N)^aB{N) (8.12) 

in complete agreement with the phenomenological Uemura relation observed in all super- 
conducting cuprate 



IX. CONCLUSIONS 

In this paper we have proposed a model of a spatially modulated collective charge state 
of the cuprates where the elementary excitations change character in real space depending 
on the local charge density. The problem is simplified by assuming that there are two main 
kind of excitations, namely, Luttinger liquid degrees of freedom in the regions of high density 
(stripes) and dx^-y'^ bound state of fermions in the regions of lower density (antiferromagnetic 
ladders), that continously transform into each other. We have shown that as a consequence 
of momentum conservation vibrations of the collective state should be present in order 
to produce sufficient phase space for condensation. In the presence of static or diagonal 
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stripes superconductivity is not possible because of the phase constraints of the former 
and the dj^2_y2 symmetry of the boson wave-function in the latter (in agreement with the 
experimental observations!!) . Therefore, unlike the BCS theory phonons are not the basic 
mechanism of pairing but without vibrations superconductivity would not be possible. 

Our main result is the phase diagram in Fig.[15|. There T* ^ Ps{N) given in ( |1.18| ) 
is the temperature scale below which dri.2-y2 bound states can form due to the confining 
potential generated by the antiferromagnetic background. This energy scale is set by the 
magnetic forces in the system and is of order of the magnetic stiffness. Because of this 
constraint our theory is only valid below the crossover line. Observe that the point where 
T* vanishes at = A^^ is not a quantum critical point as some theories assumell but it 
sets the scale for a crossover from ID to 2D behavior (see ( p. 2D ). At temperatures below 
T* boson pairs and lattice deformations start to appear in the spectrum at (±7r/a, 0) and 
(0, ±7r/a) but these deformations are dynamic in nature since long range order can only be 
attained at low temperatures due to phase fluctuations. Thus, the actual phase transition 
is driven by topological defects of the superfiuid (vortex and anti-vortex pairs) and the CCS 
distortions (dislocation loops). We have shown that at finite temperatures the transition to 
the ordered state is in the 3D XY universality class but the transition temperatures for the 
superconductivity and static lattice distortions can be very different because of the difference 
in the stiffnesses of the phase modes. In fact, we argue that static lattice distortions may be 
observable inside of the superconducting phase at very low temperatures. Although dynamic 
lattice distortions have been observed at temperatures below Tc (aseiven in ( |7.15| )) we are 
not aware of observations of static distortions at low temperatures^. The search for such 
distortions would be a good test for this theory. At zero temperature we have shown that 
phase fluctuations prevent long range order to appear until the distance between stripes reach 
the critical value A",, (given in ( |8.3| )) where phase coherence is established. For Nc < N < Ngp 
a crossover region appears in the phase diagram where incoherent bosons co-exist with stripe 
fluctuations. The phase transtion at zero temperature is in the 2D XY in a magnetic fleld 
universality class. 

So far the discussion has been based on the idea of inflnite ID stripe segments. For this 
geometry it is not possible to have a true da.2_y2 superconducting order parameter as it has 
been experimentally observecS because the Fermi surface is flat. For a flat Fermi surface 
like the one shown in Fig.^(a) the superconductivity would have s-wave symmetry since the 
whole Fermi surface would be gapped. In fact this is another major problem in approaches 
to the cuprate problem starting with ID stripes, that is, how can one get a d-wave order 
parameter out of a ID problem. In that aspect our theory indeed predicts a s-wave order 
parameter for inflnite ID stripes. If the small transverse tunneling tj_{N) given in ( ^.Ij ) is 
included, the fermion dispersion relation is modifled to ( p.3| ) and a curvature is introduced 
in the Fermi surface allowing for a very anisotropic d-wave stateil. Observe, however, that 
because t± is exponentially small with A^, changes in the Fermi surface shape will only 
become observable when A^ ^ where the system is essentially two-dimensional. This 
is not sufficient to explain the isotropic d-wave state that is observed in the underdoped 
cuprates. 

In reality, because of the presence of twin boundaries, we would expect domain formation 
with stripes running along the crystalographic directions (as it seems to be confirmed in 
neutron scattering experiments in YBC0E3). The existence of microscopic twins and tweeds 
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in the lattice structure of these systems seem to be intrinsic to the strong lattice constraints. 
This situation is quite similar to what occurs in some martensitic systemJl^. We can think 
of the twin boundaries as junctions between two superconductors. Since the twin boundary 
size, L, is very large (much larger than any of the superconducting length scales) each mono- 
domain can be treated as a separate superconductor that is coupled to other domains through 
the twin boundary. Because the coherence length in these systems is of the order of the 
lattice spacing, the superconducting order parameter can be continuously depressed near the 
twin boundary leading to a situation similar to a superconducting-normal-superconducting 
junction. As we can see in Fig.|^(a) a d^i-y^ boson couples to the positive lobe of the wave- 
function to the horizontal stripes. Thus, when phase coherence is achieved the whole Fermi 
surface in Fig.|(a) will have the same phase sign (say, positive like in Fig.P(a)). When the 
bosons cross the twin boundary and find a stripe oriented along the vertical direction the 
negative lobe of the wave-function couples to the vertical stripe. This situation is shown in 
Fig.|TB|. Therefore adjacent domains will have opposite phases across the twin boundary as 
shown in FigjT^. Thus, if we superimpose the Fermi surfaces of the horizontal and vertical 



stripes we find a situation like in Fig.|T8| where horizontal and vertical Fermi surfaces will 
have opposite sign in their phases. The situation is identitical to a vr-junction between 
two superconductor^! but the phase difference between the different domains is zero. It is 
the symmetry of the boson bound state that determines the d-wave superconductor order 
parameter. Notice that in the presence of twin boundaries the electronic motion is actually 
2D since stripe holes are being transfered from vertical to horizontal stripes via the bosons. 
Thus, there is true propagation along the diagonals. In the disordered phase this propagation 
is not coherent and the system remains essentially ID. In the superconducting state, however, 
the coherence between the domains is attained and quasiparticles can propagate along the 
diagonals as in a ordinary superconductor. Thus, our conclusion is that in the presence of 
unavoidable microscopic twin boundaries no quasi-particle peak is possible in the pseudo- 
gap phase while in the superconducting state the quasi-particle peak should exist. This 
conclusion is in agreement with the ARPES datalll. A consequence of this mechanism is 
that the size of the d-wave order parameter should change with the relative number of 
domains. In samples where there are equal amounts of vertical and horizontal domains the 



superconducting order parameter is symmetric (see Fig.|T9|(a)). If there are more domains 
in one direction than in another the d-wave order parameter should be assymetric with a 
larger lobe in one of the directions as shown in Fig.^(b). In fact, a s-wave component for 
the order parameter should also develop since for a mono-domain system we would predict 
a fully s-wave order parameter. We should also point out that c-axis coupling can lead to 
further increase of the condensation energy in the system. In fact hopping along the c-axis 
can help to stabilize the d-wave order parameteiB. 

Although Tc is weakly depedent on the twin boundaries in the sample, the transport 
properties the critical currents and fields will depend strongly on them. Let us consider a 
simple model for the junction between a vertical and a horizontal domain. The free energy 
for the junction is written as 

5F = -^cos(0) (9.1) 
where Ic is the critical current and (p is the phase difference between domains. The critical 
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current density is 



where 4 is the distance between the Cu02 planes (for simphcity we assume that Jc,a = Ic,b = 
Jc). In the presence of an electromagnetic vector potential A the Josephson current between 
two different domains is simply given by 

I J = 2e^ = h sin ((j)-^ [ dl-A] (9.3) 



$0 Jc 

where C is the line that links the two domains and $o = c/ (2e) is the magnetic flux quan- 
tum. Choosing a gauge such that = at the junction, assuming that the twin boundary 
has width d and that the variations of the order parameter across the twin boundary are 
negligible, we have from ( |9.2|) and ( |9.3| ): 

2nJrd . , 

A (9.4) 



that leads to a penetration depth, Aj, given by 



Moreover, the variation of the phase across the twin boundary is n/d but because this 
variation occurs in a length scale of the correlation length induced by the Josephson coupling, 
we expect that 

O^d. (9.6) 
Thus, the Ginzburg-Landau parameter, kj, is given by 



leading to a critical field, -ffci.j, for the field strength that is screened by surface currents 
given by 

i^.M = ^ln(..)-^ln(^-^j. (9.8) 

Finally, the critical field for the penetration of one flux of quanta through the twin boundary 
is of order 

Thus, we expect Hc2,j to be much smaller than the upper critical field required for the 
extinction of long range order in the system. It is clear that many macroscopic properties 
of the superconductor are determined by what happens at the twin boundaries. 
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As we have argued the formation of the stripes is driven by the gain in kinetic energy 
and therefore favors anti-phase domain walls as shown in Fig.^(a). Thus, the magnetic 
fluctuations are incommensurate with the lattice as it has been observed for a long time in 
all the cupratefi However, we have found that there are lattice distortions in the system 
that we can associate with the O motion as shown in Fig.^ If the stripes are site centered 
then the distortion does not affect the incommensurate spin order since the Cu atoms are 
unaffected. However, if the the stripes are bond centered (as some numerical^ as well 
as analytical work^ indicate) then the dimerization will produce commensurate magnetic 
response. This can be seen in a schematic way in Fig.^ where a dislocation loop of the CCS 
produces excess of one magnetic domain (a domain with a given staggered magnetization) 
over the otheiEl. Thus, for bond centered stripes the dislocation loops will produce a dynamic 
commensurate response. The presence of dynamical commensurate spin fluctuations have 
been observed in YBCOS but not in LSCO. It might well be that stripes are bond centered 
in YBCO and site centered in LSCO. This would explain the difference between these two 
materials in regards to the presence of commensurate spin fluctuations at low temperatures. 

It is clear from our theory that the vibrations will be affected by an isotope effect (see 
( ^.6| )). Moreover, because a static lattice distortion is involved we also expect the criti- 
cal temperature for the distortions, Tcy, to be strongly affected by changes of O isotopes 
(changes of O^^ by O^^). There is strong evidence for the O isotope effect in LSCO and other 
cuprates that support our theoryH. Moreover, we expect the isotope effect to be stronger at 
the T = transition to the superconducting state aX N = where the static lattice distor- 
tions start to appear. However, the superconducting transition itself should be very weakly 
dependent on the isotope effect because the binding mechanism only involves the lattice in 
a indirect way (in other words, the stiffness of the superfluid is only weakly renormalized 
by the lattice). The experimental data in cuprates have indeed shown an unusual isotope 
effect where the critical superconducting temperature is not correlated with isotope effect 
that becomes stronger at the quantum critical point associated with superconductivity^. 
In fact, optimally doped cuprates show very weak signs of the isotope effect. This has been 
used as an argument against phonon mechanisms for superconductivity. In our theory at 
T = a static lattice distortion appears together with the superconductivity (see Fig.^Sj). 
It is the transition into this dimerized stripe phase that will be strongly affected by the 
isotope effect, not the superconducting one. Thus, we can explain the the unusual isotope 
effect in cuprates as a consequence of the unavoidable coupling between lattice and the 
superconducting condensate. 

It is clear from the phase diagram shown in Fig.^ that our theory does not describe the 
so-called overdoped region of the cuprates. When the distance between the stripes becomes 
of the order of the lattice spacing the system becomes homogeneous. The antiferromagnetic 
correlation lenght, for instance, is short and the bound states disappear from the spectrum by 
merging with the lower Hubbard band (this coincides with the vanishing of the pseudogap 
energy scale). Beyond this point we believe there is a crossover to a conventional BCS 
behavior with a well defined Fermi surface and therefore to Fermi liquid behavior in the 
normal phased!. This smooth crossover is possibly the same one that occurs between a 
Bose-Einstein system and a BCS superconductorl^. 

In summary, we have presented a new model for a collective electronic state of the 
cuprates where the elementary excitations change from place to place in real space. We show 
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that the decaying processes among these elementary excitations produce superconducting 
correlations even when the interactions are repulsive. We show that the d-wave nature of 
the order parameter is associated with the d-wave nature of boson bound states that exist 
due to the magnetic confinement in the system. We have shown that phase fluctuations are 
responsible for the quantum disorder and that the phase diagram depends strongly on how 
vortices couple to dislocation loops. We have explained various different experimental facts 
of the cuprates and predicted new effects that might prove or disprove our theory. 
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FIG. 1. Schematic plot of the density of states as a function of energy. 




FIG. 2. Effective potential felt by the carriers: (a) K(y) is the string potential generated by 
the antiferromagnetic background and Va(2/) is the atomic potential generated by the lattice; (b) 
resultant potential felt by the carriers. 
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FIG. 3. (a) Magnetic interactions in the presence of stripes: J is the exchange between spins 
in the antiferromagnetic ladder; J' is the exchange across the stripe, (b) Effective magnetic model 
with spatially anisotropic couplings. 
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FIG. 4. Momentum space picture for stripes: (a) longitudinal; (b) diagonal. 
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FIG. 5. Basic decaying processes for static stripes: (a) momentum exchange; (b) real space 
process where a boson pair can decay into stripe fermions. 
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FIG. 6. CoupUng of a dr^2_y2 boson wavefunction to (a) longitudinal and (b) diagonal stripes. 
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FIG. 7. Scattering processes including vibrations. 



[•"o'iT o •"o"ir o •"o"ir o "(ij 

o'"'o""" o"""o""" o"""o""" 

•o«o«o*o«o«o* 

0000000 

•o«o«o«o«o«o* 

0000000 

•o»o«o«o»o*o« 

0000000 
!«'6'» o •"o"« o •"o"» O 

o"""o""" o"""o""" o"""6""" 

•0*0*0«0«0*0« 

0000000 
•o«o*o«o«o«o« 

0000000 
•o«o«o«o«o«o« 

0000000 

FIG. 8. Real space picture of the normal phase. Open circles: O atoms; Filled circles: Cu 
atoms; dashed line: equilibrium position of the stripes. 
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FIG. 9. Real space picture of the ordered phase: 5u indicates the possible motion of the O 
atoms with doubling of the unit cell perpendicular to the stripe orientation. The difference in gray 
levels represents the modulation of the superfluid density. 
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FIG. 10. Momentum space picture of the condensation region. 
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FIG. 13. Critical value of the coupling constant F as a function of 1/iV. 
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FIG. 14. Topological defect of the CCS. 
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FIG. 15. Phase diagram of the model (the symbols are explained in the text). 
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FIG. 16. Coupling generated by a dj.2_y2 boson in the presence of vertical and horizontal stripes. 
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FIG. 19. (a) Superconducting order parameter for equal amounts of vertical and horizontal 
domains; (b) superconducting order parameter in the case where there is an excess of domains in 
one direction. 
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FIG. 20. + si 
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netization: (a) normal phase; (b) dimerized phase; (c) topological defect of a bond centered stripe. 
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